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Periodicity is one of the most fundamental structural characteristics of systems occurring in 
nature. The properties of these systems depend strongly on the symmetry of the underlying periodic 
structure. In solid state materials - for example - the static and transport properties as well as 
the magnetic and electronic characteristics are crucially influenced by the crystal symmetry. In 
this context, hexagonal structures play an extremely important role and lead to novel physics like 
that of carbon nanotubes or graphene. Here we report on the first realization of ultracold atoms 
in a spin-dependent optical lattice with hexagonal symmetry. We show how combined effects of 
the lattice and interactions between atoms lead to a forced antiferromagnetic Neel order when two 
spin-components localize at different lattice sites. We also demonstrate that the coexistence of two 
components - one Mott-insulating and the other one superfluid - leads to the formation of a forced 
supersolid. Our observations are consistent with theoretical predictions using Gutz wilier mean- field 
theory. 
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In recent years, ultracold atoms in periodic 
potentials^ ^ have been widely recognized as an im- 
portant tool to simulate solid state systems and study 
their transport^ and magnetic properties^ ^ . Different 
types of spin- and/or state-dependent lattices have 
been implemented and studiecP^^. So far, most exper- 
iments with ultracold quantum gases have been carried 
out in lattices of cubic ^mmetry. However, recent 
theoretical development^^^"^ are aiming at systems 
with a hexagonal geometry. In particular, carbon 
nanotube^I^, graphene^^ and a large number of other 
carbon-based compounds show fascinating effects and 
exhibit particularly rich quantum phases^. In this paper 
we discuss the first realization of ultracold quantum 
gases in a hexagonal, spin-dependent optical potential 
and demonstrate how it can be used to tailor quantum 
phases of spin-mixtures and their dynamics. We show 
that the combination of interactions between different 
spin-states and the spin-dependent lattice potential 
leads to novel quantum phases: a forced antiferromag- 
netic Neel order when two spin-components localize in 
different sublattices, and a forced supersolid phase, with 
one spin-component being localized in one sublattice and 
inducing a density modulation on the other superfluid 
spin-component. These phenomena are studied by 
exploiting a novel technique of state- and site-selective 
microwave-spectroscopy. Examining the impact on the 
onset of the superfluid to Mott-insulator transition of 
the respective spin- state, we demonstrate furthermore 
that the mobility of particles in the lattice can be 
adjusted by immersing a well localized spin-crystal into 
a superfluid bath. We calculate the corresponding phase 
diagrams within Gutzwiller mean-field theory and find 
good agreement with the experimental results. 



I. GENERATION OF A HEXAGONAL, 
SPIN-DEPENDENT OPTICAL LATTICE 

The basic structure of the spin-dependent hexagonal 
optical lattice discussed in this paper is illustrated in 
Fig. [T] Three laser-beams intersecting under an angle of 
120° and with each beam linearly polarized in the plane 
of intersection lead to the formation of local potential 
minima in a hexagonal structure. For neighboring sites 
along the vertices of the hexagonal lattice as for example 
indicated in Fig. [T^ by a dashed line, the local polariza- 
tion alternates between cr^ and (j~ . As atoms in a light 
field experience a polarization-dependent ac Stark shift, 
the potential at these <j+ and <j~ sites is different for 
different atomic Zeeman substates labeled by mF- The 
potential can be written as: 



l/(x) = 14ex (x) + rriFgFl^BBeE (x) , 



(1) 



where the polarization of the light field V (x) (+1 for pure 
(7+ and —1 for pure (j~ polarizations) is mapped onto 
a pseudo-magnetic field 5eff (x) oc — Vhex (x) V (x) //ie, 
qf the Lande ^-factor and /iB the Bohr magneton (for 
details see appendix). The potential consists of a spin- 
independent part Fhex (x) of hexagonal symmetry and 
a state-dependent super-lattice emerging from the local 
pseudo-magnetic field B^f^ (x). According to the local po- 
larization, we denote the emerging triangular substruc- 
tures as (7+ and g~ lattice (see Fig. Ip). The hexagonal 
lattice can therefore also be regarded as a triangular lat- 
tice with a bi-atomic basis where the atoms occupy a^ 
and G~ sites as indicated in Fig.jlp by green and red bul- 
lets. For typical experimental parameters, the mobility 
of the atoms is predominantly restricted to the hexag- 
onal channel structure connecting the local minima of 
"^hex (x) . One possible tunneling channel is indicated in 
Fig. [T^ by a dashed line. Note that the large central 
peak maxima (9l^at) prevent a direct diagonal tunneling 
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Figure 1. a Hexagonal optical lattice with alternating a~^ (green bullets) and a~ (red bullets) polarization structure. The 
upper ID plot shows a cut through the 2D potential. The tunneling of atoms loaded into the potential is mainly restricted 
to channels like the one indicated by the orange dashed line. The lowest graph shows the 2D polarization distribution in 
the lattice ranging from fully cr^ (green) to cr~ (red) polarized. The hexagonal structure is based on two triangular sublattices 
which are translated with respect to each other, b One-dimensional potential structure along the channel shown in subfigure 
a for particles in different Zeeman states |F, ttif). The modulation depth f{mF)Viat of the channel depends on the m^-state 
and is given by {l,1.7,2.5}yiat for ttif — {0,=bl,zb2} and A = 830nm. The degeneracy of neighboring lattice sites for atomic 
substates with ttif / is lifted by the polarization- and atomic state-dependence of the ac Stark shift. The lower part shows 
the corresponding light polarization. When the potential is deep enough, atoms in the different substates will localize as 
indicated, c Time-of- flight images of superfluid samples prepared in different Zeeman substates released from the lattice (left) 
and corresponding calculated quasi- momentum distributions (right). The second- to first- order interference peak ratio changes 
from the triangular to the hexagonal lattice occupation which is in direct agreement with the different structure factors of the 
two lattice geometries. 



through these maxima. Figure [TJd illustrates the different 
lattice potentials along a tunneling channel for the F = 1 
hyperfine ground-state manifold of ^^Rb. The F = 2 hy- 
perfine states exhibit a similar behavior with twice the 
modulation amplitude for \mF\ = 2 (compare Eq. ([l])). 
To study atomic quantum phases in the lattice, we load 
an ensemble of ultracold ^^Rb ground-state atoms into 
the lattice. This consists of either pure or specific com- 
positions of different hyperfine- and magnetic Zeeman- 
states. Due to the additional spin-dependent spatial vari- 
ation of the potential, atoms with tuf y^ are preferably 
confined to either the <j+ or the a~ triangular lattice 
structure. In contrast, atoms in the tuf = state will 
distribute homogeneously over both the a~^ and a~ lat- 
tices (Fig. lib). For convenience, we refer to these differ- 
ent configurations as triangular {tuf y^ 0) and hexagonal 
states {rriF = 0). This spin-dependent spatial distribu- 



tion determines the structure factor of the lattice, which 
in turn determines the quasi-momentum distribution of 
the atoms. We experimentally observe this difference in 
the structure factors by mapping the quasi-momentum 
spectrum to the spatial density distribution of the ensem- 
ble via time-of-flight images. Typical results for triangu- 
lar and hexagonal initial states are shown in Fig. [T]: to- 
gether with the corresponding theoretical prediction from 
a one-particle band structure calculation. The observed 
second- to first-order interference peak ratio (see Fig. [l]: 
bottom) is significantly larger for the hexagonal states 
as compared to triangular configurations. It typically 
ranges between two and four for our experimental pa- 
rameters, which is well reproduced by calculations of the 
Bloch functions and their quasi-momentum distribution 
for the corresponding spin-dependent lattice geometry. 
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Figure 2. Results of microwave-spectroscopy measurements. The addressed transitions are indicated in the insets. For com- 
parison, the lattice depth V^iat is given for each data set. Graphs a-c,e show measurements for single-component samples and 
graph d for a spin-mixture loaded into the lattice. The frequency axis is shifted by the resonance frequency cjzero of the re- 
spective sample for vanishing 2D-lattice (Mat = 0). a Transferring the atoms from a hexagonal state |1,0) to another hexagonal 
state |2,0) leads to a single Fourier limited transition line at cjq- b Spectroscopy of the triangular state |1, —1) to |2,0). For 
sufficiently deep lattice potentials, the atoms start to localize in the deeper wells. The transition spectrum therefore shows a 
dominant resonance at the corresponding frequency cji. c Addressing the |1, 0) -^ |2, — 1) transition leads to two well separated 
resonance lines UJ2 and cjs corresponding to the energy splitting of the potential wells due to the pseudo-magnetic field Beff. 
The integrals of the individual peaks are proportional to the number of atoms in the respective wells. The different height of 
the peaks reflects the different Franck- Condon overlap between the initial and the different final states, d Interaction induced 
localization in spin-mixtures. An equal mixture of |1,0) and |1, — 1) atoms is prepared and the occupation in the lattice wells 
is determined for the |1,0) atoms. The presence of |1,— 1) atoms suppresses the population of the |1,0) atoms at a~ sites 
and only the a~^ resonance remains, e Spectroscopy data for different lattice depths. The peak centered at corresponds to 
the transition frequency in the system without lattice. For increasing lattice depth, the atoms redistribute and localize in the 
deeper wells. 



II. STATE- AND SITE-SELECTIVE 
MICROWAVE SPECTROSCOPY 



To obtain a detailed understanding of the rich physics 
arising from different lattice occupations we have devel- 
oped a novel method for state- and site-selective mi- 
crowave spectroscopy. It allows for in-situ investiga- 
tions of the spatial spin-ordering in the hexagonal and 
triangular lattice structures. It is suitable for study- 
ing different regimes ranging from the superfluid to 
the Mott-insulating phase. Previously, a different kind 
of polarization-dependent radio-frequency spectroscopy 
method was used to study cubic lattices^. Our method 
is based on the spatially varying transition frequencies 
between populated 772i?-states to initially non-occupied 
TTiF' -states. The potential energy difference between 



the initial and the final state is given by AV(x) = 
Mb {rriF' + tuf) {Bq + 5eff (x)) /2 + ^^hfs, where Bq is an 
additionally applied homogeneous magnetic guiding field 
and ^hfs is the ^^Rb- hyper fine splitting. The spec- 
troscopy signal is obtained by detecting and normalizing 
the number of atoms in the respective rriF- and tuf'- 
states after applying a microwave pulse, then releasing 
the sample from the lattice and performing Stern- Gerlach 
separation. In the following we discuss central results of 
the microwave spectroscopy method, first for single com- 
ponent samples, then for spin-mixtures. Typical data for 
both scenarios is shown in Fig. |2] 

The spatial variation of the differential potential shift 
AV(x) vanishes for transitions with tuf' -^ttif = 0. Thus 
a single component sample initialized in the \F^mF) = 



1 1,0) state shows a single transition line when the state 
1 2,0) is addressed (Fig. [2^). When transferring atoms 
from the hexagonal state"l, 0) to a triangular final state 
{rriF' -\-raF 7^ 0; the data in Fig. [2]: shows |1, 0) to |2, — 1) 
for Flat = S.TE^rec), the spectrum splits up into two lines 
according to the energy difference between the neighbor- 
ing lattice sites. Here, the left (right) transition peak 
refers to the transition from the first band of the hexago- 
nal stdite |1, 0) to the first (second) band of the triangular 
state |2, — 1) where the atoms are mostly confined to the 
a~{a~^) sites. We measure the two transition frequencies 

6^2^^ — CJzero = 27r X ( — 3.6 ± 0.2)kHz and CJg^^ — CJzero = 

27rx(7.1±0.2)kHz, where cjzero is the resonance frequency 
without 2D-lattice. We find excellent agreement with the 
theoretical expectations of co'2^®° — ^zero = — 27rx3.8kHz 
and co'3^®° — cc;zero = 27rx7.1 kHz obtained using an ab ini- 
tio single-particle 2D band-structure calculation. Note 
that for the hexagonal lattice the Bloch functions do not 
separate like it is the case for cubic lattice. 

Moreover, we compare the particle transfer efficiencies 
of the two transitions with the calculated predictions. 
Experimentally, the transfer efficiencies are determined 
evaluating the integrals of the corresponding transition 
peaks. For the measurement shown in Fig. [2]: we measure 
the transfer efficiency of the a;2-transition to be (2.0±0.2) 
times larger than for the cJs-transition. This is in very 
good agreement with the calculated ratio of the Franck- 
Condon overlaps between the initial state and the two fi- 
nal states resulting in a 1.9 times larger transfer efficiency 
for the CJ2 -transition. It also confirms the expected ho- 
mogeneous density distribution of tuf = atoms at a~^ 
and (j~ sites. Therefore the spectroscopy allows to mea- 
sure the relative occupations of the two sublattices cr+ 
and a~ . It can equally well be applied to study other 
initial spin-states e.g. a state of |1,— 1) atoms. Figures 
^p and e show the |1, — 1) to |2,0) transition amplitudes 
for different lattice depths and clearly reveal the exclusive 
occupation of the a~ sites for sufficiently deep lattices. 
For shallow lattices, the microwave transition frequency 
coincides with the one for a vanishing lattice potential. A 
crossover can be observed at intermediate lattice depths. 
The observed shift of 5 — 7kHz for deeper lattices is repro- 
duced by calculations within a harmonic approximation, 
whereas the signals for shallow and intermediate lattices 
cannot easily be explained in the single-site approxima- 
tion. When initializing the atoms in |1,+1), we observe 
the equivalent signals corresponding to the occupation of 
the (7+ sites and thus demonstrate the ability to actively 
force the generation of a magnetically ordered Neel state 
with predefined magnetization by using a mixture of dif- 
ferent spin-states, e.g., |1,— 1) with |1,+1). 

The experimental setup allows for the preparation of ar- 
bitrary spin-mixtures in the spin-dependent optical lat- 
tice. The interplay of the interaction between different 
spin-states and the lattice potential has a crucial impact 
on the static and dynamic properties of the system. Par- 



ticularly interesting is the case of a mixture of |1, 0) and 
|1,— 1) atoms. While the |1,— 1) atoms tend to localize 
in the a~ sites, the |1,0) atoms fill the hexagonal lattice 
homogeneously as it was shown above for pure atomic 
samples (Figs. [2^-c, e). In case of a mixture of these 
spin-components, the additional repulsive m^er-species 
interaction causes the |1, — 1) atoms to imprint a periodic 
density modulation on the |1,0) atoms. Experimentally 
we observe an almost vanishing transition amplitude of 
the peak corresponding to the a~ sites as shown in Fig- 
ure [2]i. Comparison of the integrated overall transition 
probability with the one obtained for pure systems in- 
dicate that (30 ±5)% of the tuf = atoms, initially 
sitting at a~ sites, have been transferred to a~^ sites. 
This crystalline order imprinted onto the |1,0) atoms by 
the 1 1,-1) atoms can appear in two different scenarios: 
First, when the |1, 0) atoms are localized on the (T+ sites, 
the system forms an alternating Mott-insulator. In this 
case, the mixture resembles a localized alternating spin- 
ordering. Second, when the |1,0) atoms are still super- 
ffuid the system shows an off-diagonal long-range order. 
Together with the spatial modulation due to the interac- 
tion with the other s pin-co mponent , this indicates a su- 
persolid like behavioi^^^EH. By varying the depth of the 
lattice both scenarios have been realized and indeed, we 
experimentally observe that the |1,0) atoms were driven 
across the superffuid-to-Mott-insulator transition refer- 
ring to both cases. 

Our observations are strongly supported by a theoretical 
analysis using the mean-field Gutzwiller approximation^ 
(see Appendix). In Fig. Islwe demonstrate that an equal 
mixture of |1, 0) and |1, — 1) atoms results in a forced su- 
persolid. This is shown with the aid of the absolute value 
of the superfiuid order parameter (a^ri) (top) and the 
occupation number (n^ri) (bottom) of the |1,0) state at 
<j+ and (j~ sites assuming a filling of two particles per 
unit cell. For comparison, we also include the calcula- 
tions for a pure |1,0) system. Here, | {aa+) I and | (a^r-) | 
as well as {h^+) and {n^r-) coincide for all Fiat reveal- 
ing the expected homogeneous atom distribution within 
the lattice. This changes when |1,— 1) atoms are ad- 
mixed to the 1 1,0) atoms. At already moderate lattice 
depths the |1, — 1) particles are localized at the cf~ sites. 
As a consequence the |1,0) component is repelled from 
these sites and the corresponding occupation number and 
superfiuid order parameter become spatially modulated, 
which is characteristic for a supersolid. Moreover, we 
observe that the effective blocking of (j~ sites by |1, —1) 
particles forces the |1, 0) atoms to undergo the transition 
to the Mott-insulator already at weaker lattice depths 
than in the pure case. Hence an admixture of a second 
spin-component allows us to tune the critical point of the 
SF-MI transition as well as the spatial density modula- 
tion and superffuid order parameter. 
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Figure 3. Results of Gutzwiller mean- field theory for a) occu- 
pation number (n^±) and b) absolute value of the superfluid 
order parameter | (a^±) | of |1,0) atoms. Two different en- 
sembles are considered; a pure ensemble of |1,0) atoms (blue 
line) and a 50:50 mixture with additional |1, — 1) atoms (black 
lines). The solid (dashed) black line is the value at g^ [^'^~^ 
sites. In a pure system the values are identical at g^ and g~ 
sites (blue line). The particle number of each spin-state is 
fixed to two per unit cell. 



III. SUPERFLUID TO MOTT-INSULATOR 

TRANSITION IN A SPIN-DEPENDENT 

OPTICAL LATTICE 



Transport properties play an important role in under- 
standing, for example, the conductivity of solid-state sys- 
tems. For the quantum optical counterparts of such sys- 
tems, the optical lattices, these properties are essentially 
governed by the interplay of the on-site interaction U 
with the the tunnel matrix element J and are usually well 
described by the Bose-Hubbard modeP (see Appendix) . 
In optical lattice experiments, the ratio of J jU is di- 
rectly adjustable through the power of the lattice beams 
allowing to drive the ensemble from the superfluid to the 
Mott-insulating state. We use the interference pattern 
contrast (visibility) as an indicator for this transition^^. 
Results of such measurements are shown in Figure [4^ 
for different single-component samples. We compare the 
SF-MI transition of a hexagonal state |1, 0) with that of 
the triangular states |1, — 1) and |2, — 2) and observe that 
the interference contrast curves are considerably shifted 
with respect to one another. 

We start by analyzing the situation for the |1, — 1) and the 
1 2, —2) states. As these states experience an additional 



modulation of the hexagonal potential caused by B^f^ (x) , 
resonant tunneling is only possible to the six next-nearest 
neighbor sites in the lattice. Due to this additional modu- 
lation, the tunneling barrier between the occupied lattice 
sites increases with Im^l (compare Eq. Q). The SF-MI 
transition is thus expected to occur at lower overall beam 
powers for increasing values of |mi?|. However, the two 
transitions are expected to nearly coincide when both 
are rescaled to the same modulation depth. This is con- 
firmed by the data shown in the inset of Figure [4^, where 
the interference contrast is plotted as a function of the 
effective one-dimensional tunneling barrier height. For 
the rriF = states, the situation is different: Since in 
this case the optical potential is not additionally modu- 
lated, resonant tunneling takes place to the three nearest 
neighbors. This leads to a strongly increased tunneling 
probability and therefore higher lattice beam powers are 
required for the SF-MI transition to occur. Interestingly, 
all three transitions nearly coincide when normalized to 
the height of the barrier, regardless of the effect of differ- 
ent hopping distances and possible number of accessible 
sites. This discrepancy can be explained when consid- 
ering in addition to nearest neighbor tunneling processes 
also second-nearest neighbor processes. This is confirmed 
in calculations using the Gutzwiller approach. Results of 
such calculations are shown in Fig. ^p for the experimen- 
tal configurations with filling factors ranging between two 
and four per lattice unit cell. The numerical results show 
a good qualitative agreement with the observed experi- 
mental data. 

The fact that the SF-MI transition occurs at different lat- 
tice beam intensities for different values of \mF\ adds an 
exciting new dimension to our lattice geometry: it opens 
the possibility to create spin-mixtures where a fully local- 
ized component is immersed in a superfluid bath. Exper- 
imentally, we study this configuration by loading |1, — 1) 
atoms which preferably occupy the a~ lattice and add to 
these the more localized state |2, — 2) (|2,+2)), which is 
strongly confined to the a~^ {(^~) lattice. As Figure Ep 
shows for the case of a 50:50 mixture, a considerable shift 
in the SF-MI transition of the |1, — 1) atoms is observed 
when the admixed atoms occupy a~^ sites, i.e., when they 
are located in between the preferably occupied a" sites of 
|1, — 1). On the other hand, the transition remains prac- 
tically unaltered when the more localized component oc- 
cupies the same lattice sites as the superfluid core. The 
pictograms in Figure [4J3 illustrate these three different 
cases. The addition of a spin-crystal thus allows one to 
tailor the position of the SF-MI transition depending on 
the sub-lattice it occupies. 

The Gutzwiller computations reproduce the shift of the 
SF-MI transition (see Figure |4]i). The observed changes 
in the tunneling probabilities can be ascribed to a modi- 
fication of second order hopping processes: Second order 
hopping of |1,— 1) between the <j~ sites is only possi- 
ble by passing through sites of the a~^ lattice. However, 
in a |1,— 1) - |2, — 2) mixture, these sites are occupied 
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Figure 4. Experimental observation of tlie superfluid to Mott-insulator (SF-MI) transition for pure and mixed atomic spin-state 
system averaged over typically 7-10 experimental runs (a and c) and the corresponding theoretical superfluid order parameter 
normalized to the square root of the total occupation number for the same configurations, calculated within a mean- field 
Gutzwiller approach (b and d). a Visibility curves for different single spin-components of the optical lattice as function of the 
lattice depth Viat, which is proportional to the lattice laser beam powers. In the inset, the data is presented with the lattice 
depth renormalized to the height of the nearest neighbor tunneling barrier, b Theoretical calculation of the superfluid order 
parameter normalized to the square root of the total occupation number for the experimental settings presented in a. The 
approach of the superfluid order parameter to zero indicates the position of the SF-MI transition. The areas represent the 
results for lattice filling factors between two and four per unit cell as expected for our experimental parameters, c Comparison 
of visibility curves for pure |1, —1) and 50:50 mixtures of |1, —1) with |2, —2) (|2, +2)) atoms. The insets indicate the preferred 
probability distributions of the atoms for the different mixtures and the modification of the tunneling matrix element J of 
|1, — 1). d Theoretical calculation of the superfluid order parameter of the |1, — 1) state for filling factors between two and four 
per unit cell. 



by |2, — 2) atoms. The interaction energy between the 
atoms increases the effective barrier and suppresses tun- 
neling. Thus, the SF-MI transition is shifted towards 
lower lattice laser beam intensities. By changing the rel- 
ative composition of the mixture, the strength of this 
tunneling suppression can be modified. 
Further insight into the underlying physical processes is 
gained by considering the phase diagrams calculated in 
Gutzwiller approximation. Figure [5^ shows the SF-MI 
phase diagram for the experimental configurations dis- 
cussed in Figs. [4J: and d. A clear modification of the 
Mott-insulator lobes of a triangular state can be observed 
when strongly localized atoms are present in the comple- 
mentary triangular sublattice(|l, —1) + |2,— 2)). Such a 
strong modification with respect to the pure case does 
not occur when both components occupy the same sites 
(|1,— 1) + |2,+2)). This perfectly agrees with our ex- 
perimental observations discussed above. In Fig. l5b we 



plot the SF-MI transition for a |1,-1) - |2, -2) mix- 
ture for different inter-component interaction strengths 
U-i-2 between |1,— 1) and |2, — 2) ranging from the 
experimental value /7_i^_2 to 30/7_i^_2. Such a sit- 
uation of e.g. 30/7-1,-2 could be achieved experimen- 
tally by tuning the interactions with the help of Fesh- 
bach resonances^ ^. Interestingly, the standard tendency 
is inverted in the hard core limit, when different compo- 
nents occupy different sites. This implies that the SF-MI 
transition for higher Mott orders occurs in the interest- 
ing regime that interpolates between the Bose-Hubbard 
regime and an array o f coup led rotors, similar to an array 
of Josephson junctions 
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Figure 5. Superfluid to Mott-insulator transition phase diagram from Gutzwiller mean-field theory for |1,— 1) atoms as a 
function of J-i/V^i and the chemical potential /i_i/[/^i . a An admixture of |2,— 2)-atoms (blue dot-dashed line) strongly 
enlarges the |1,— 1) Mott-lobes with respect to the pure system (red solid line) while an admixture of |2,+2)-atoms (green 
dashed line) hardly affects the transition, b In a |1, — 1)-|2, — 2)-mixture the |1,— 1) Mott-lobes grow when increasing the 
interaction strength between the spin-components (given in multiples of the ^^Rb atomic interspecies interaction strength 
between |1, -1) and |2, -2)). 



IV. CONCLUSION 
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In conclusion, we have presented a comprehensive ex- 
perimental and theoretical study of a novel type of optical 
lattice, which exhibits hexagonal spatial symmetry. The 
intriguing polarization structure of the lattice sites can 
cause a break up of the hexagonal structure for partic- 
ular spin-states leading to the creation of a triangular, 
magnetic sublattice. The nature of this lattice potential, 
combined with the possibility to precisely control sys- 
tem parameters, makes this system especially suitable 
for studying new quantum phase mixtures of superfluid 
and Mott-insulating states. 

We could show that a combined effect of repulsive interac- 
tion and the specific hexagonal lattice structure leads to 
modifications of J//7, which is clearly a different situation 
compared to e.g. Bose-Fermi mixture experimenta^^^ 
where higher order effects effectively modify J and U . 
As promising further directions, entropy effects and en- 
tropic cooling^^ can be studied in the hexagonal sys- 
tem. Particularly interesting, however, is the prospect of 
studying polarized or unpolarized Fermi gases, or Fermi- 
Bose mixtures in the hexagonal lattice geometry. At half 
filling, both of these polarized systems will mimic the 
physics of graphene^, i.e., exhibit Dirac's dispersion re- 
lation, artificial relativistic effects, etc. Introducing ar- 
tificial gauge fields in such situations^ ^^ will lead to the 
occurrence of the anomalous quantum Hall effect and a 
whole variety of exotic quantum phase transitions (see 
for instance Ref.^*^). Furthermore, by only changing the 
polarization of the laser beams creating the lattice, a 
spin-independent triangular lattice can be created^. A 
bosonic gas in the triangular lattice might then be used 
to mimic frustrated antiferromagnetism by employing a 
time-dependent lattice modulation^^. 
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APPENDIX 

Spin-dependent hexagonal optical lattice. The 

light for the optical lattice is derived from a Ti: sapphire 
laser running at A = 830 nm. To realize the two- 
dimensional hexagonal lattice we interfere three laser- 
beams under an angle of 120° in one plane. The po- 
larizations of the beams are linear and in the plane of 
the intersecting area and perpendicular to a homoge- 
neous magnetic guiding field of Bq = I.IG which de- 
fines the quantization- axis. In order to compensate for 
phase- noise, the phase of each laser-beam is actively sta- 
bilized. To ensure two-dimensionality of the hexagonal 
structure, a perpendicular one-dimensional lattice with 
different frequency is applied. Typical potential depths 
of this ID-lattice are ^id = 44^rec- The resulting two- 
dimensional, hexagonal lattice potential is given by 

V(X) = -2Fiat{3 - COS [(bi - bs) • X] + COS (bi • x) 

+ COS (b2 • x) I - 2V3Vis,t0^mF\ sin [(bi - b2) • x] 
+ sin (bi • x) — sin (b2 • x) > 

= Fhex (x) + mF^FMB^eff (x) , (2) 

where bi = ^/3/2kex + k/2ey and b2 = V^ke^ with k = 
27r/A are the reciprocal vectors of the hexagonal lattice. 
The constant a depends on the detuning of the lattice 
laser from the ^^Rb resonances at 780/795nm and is given 
for the experiments reported here (A = 830 nm) a = 0.13. 
The total potential is then given by V^(x) + Vid sin^(/cz). 
The one- dimensional representation (Figure III)) can be 
reduced to a sinusoidal potential with 319 nm period, 
modulated by a cosine function with 638 nm period. The 
energy structure of the system can be readily calculated 
using an ab initio single-particle 2D band-structure cal- 
culation. 

Preparation scheme. We start with a Bose-Einstein 
condensate of typically several 10^ atoms in the stretched 
state 1 1,-1), which are confined in a nearly isotropic 
crossed dipole trap (cj ~ 27r x QOHz). The preparation 
of the different pure and mixed spin-states is performed 
with aid of radio-frequency and/or microwave sweeps. 
After the state preparation we ramp up the optical lat- 
tice within 80ms using an exponential ramp. Within the 
ramping time the coherence between different spin-states 
is lost. 

Spectroscopy. We apply a microwave square-pulse of 
typically 2 ms duration allowing for a sufficient resolu- 
tion in frequency space. In order to minimize finite- 
state interactions, the intensity of the microwave pulse 
is adjusted such that typically only a small fraction of 



up to 20% of the atoms populate the final state. To 
separate different spin-components during 27 ms time-of- 
flight prior to absorption imaging, a Stern-Gerlach gra- 
dient field is applied. 



Bose— Hubbard model and Gutzwiller mean-field 
Ansatz. For the computation of the theoretical phase 
diagrams (Figs. [IJd, d, and Fig. |5| the system is modeled 
by a Bose-Hubbard model, which is known to give an es- 
sentially exact description of ultra-cold atoms in optical 
lattices^. The corresponding Hamiltonian reads 



H= ^ { - J, ^ (d|d,- + h.c.) 

a=a,b (ij) 

+ 2 ^ (^" K'^ ("""'^ " ^^ " ('"" " ^«)^"'^ 



(3) 



where a and b denote the two different spin-components 
and a~^ and a~ distinguish between the two sublattices. 
The operator a^ (bi) destroys an a (b) boson at site i, and 
^a,i {^h,i) is the corresponding occupation operator. An- 
gle brackets denote pairs of nearest neighbors. The tunnel 
matrix element is denoted by J^^ the on-site interaction 
between the spin-components \a) and \a') at a^ site by 

^a,t'(^a^=^a,t)- Local potentials are denoted by E^ 
and the chemical potentials by jIq^. The parameters £"2, 
/7^Q,/, and Jo^ are extracted from the experimental set- 
ting. For the derivation of the on-site interaction U we 
integrate over the Gaussians which approximate well the 
Wannier functions at the minima of the optical lattice. 
We estimate the nearest neighbor tunneling J from exact 
one-particle band-structure calculations of the full optical 
lattice. Comparing this to a tight-binding Hamiltonian 
gives the desired tunneling coefficients. 
To solve the model Hamiltonian Eq. (Isl) we have em- 
ployed the standard Gutzwiller mean-field approxima- 
tion which is known to give good qualitative and quan- 
titative account of the phase diagram in two and more 
dimensions and which reproduces the exact solution for 
U/J ^ oo and, in the thermodynamic limit, also for 
U/J -^ 0.^ The superfluid order parameter is defined as 
the expectation value of the annihilation operator. In 
order to calculate the phase diagrams in Fig. [5J we have 
chosen the chemical potential to be equal for both com- 
ponents, i.e., /i_i = /i±2- In the |1, —1) — |2, —2) mixture 
this leads to an approximately 50:50 mixture, while in 
the |1, —1) — |2, +2) mixture the |1, —1) atoms repel the 
1 2, +2) strongly such that the occupation number can be 
heavily asymmetric. 



sengstock@physik.uni-hamburg.de' 
^ Jaksch, D., Bruder, C, Cirac, J. I., Gardiner, C. W., and 

ZoUer, P. Phys. Rev. Lett SI, 3108 (1998). 
^ Greiner, M., Mandel, O., EssUnger, T., Hansch, T. W., 

and Bloch, I. Nature 415, 39 (2002). 
^ Trotzky, S., Cheinet, P., Folhng, S., Feld, M., Schnor- 

rberger, U., Rey, A. M., Polkovnikov, A., Demler, E. A., 

Lukin, M. D., and Bloch, I. Science 319(5861), 295-299 

(2008). 
^ Lewenstein, M., Sanpera, A., Ahufinger, V., Damski, B., 

Sen, A., and Sen, U. Adv. Phys. 56, 243 (2007). 
^ Bloch, I., Dalibard, J., and Zwerger, W. Rev. Mod. Phys. 

80, 885 (2008). 
^ Jaksch, D., Briegel, H. J., Cirac, J. I., Gardiner, C. W., 

and Zoller, P. Phys. Rev. Lett 82, 1975 (1999). 
^ Mandel, O., Greiner, M., Widera, A., Rom, T., Hansch, 

T. W., and Bloch, I. Nature 425, 937 (2003). 
^ Lee, P. J., Anderlini, M., Brown, B. L., Sebb-Strabley, J., 

Phillips, W. D., and Porto, J. V. Phys. Rev. Lett. 99, 

020402 (2007). 
^ Lin, Y.-J., Compton, R. L., Perry, A. R., Phillips, W. D., 

Porto, J. v., and Spielman, L B. Phys. Rev. Lett. 102, 

130401 (2009). 
^° Lin, Y.-J., Compton, R. L., Jimenez-Garcia, K., Porto, 

J. v., and Spielman, L B. Nature 462, 628 (2009). 
^^ McKay, D. and DeMarco, B. arXiv:0911.4143v, to appear 

in New J. Phys. (2010). 
^^ Eckardt, A., Hauke, P., Soltan-Panahi, P., Becker, C, Sen- 

gstock, K., and Lewenstein, M. Europhys. Lett. 89, 10010 

(2010). 
^^ Snoek, M. and Hofstetter, W. Phys. Rev. A 76(5), 051603 

Nov (2007). 



^^ Lee, K. L., Gremaud, B., Han, R., Englert, B.-G., and 

Miniatura, C. Phys. Rev. A 80(4), 043411 (2009). 
^^ Geim, A. and Novoselov, K. Nature Materials (2007). 
^^ lijima, S. Nature 354, 56 (1991). 
^^ Du, X., Skachko, L, Duerr, P., Luican, A., and Andrei, 

E. Y. Nature 462, 192 (2009). 
^^ Biichler, H. P. and Blatter, G. Phys. Rev. Lett. 91(13), 

130404 (2003). 
^^ Titvinidze, L, Snoek, M., and Hofstetter, W. Phys. Rev. 

Lett. 100(10), 100401 (2008). 
^° Gerbier, P., Widera, A., Foiling, S., Mandel, O., Gericke, 

T., and Bloch, L Phys. Rev. A 72(5), 053606 (2005). 
^^ Inouye, S., Andrews, M. R., Stenger, J., Miesner, H.-J., 

Stamper-Kurn, D. M., and Ketterle, W. A^a^i^re 392(6672), 

151-154 (1998). 
22 Doniach, S. Phys. Rev. 5 24(9), 5063-5070 (1981). 
2^ Sachdev, S. Quantum Phase Transitions. Cambridge Uni- 
versity Press, Cambridge, (1999). 
2^^ Giinter, K., Stoferle, T., Moritz, H., Kohl, M., and 

Esslinger, T. Phys. Rev. Lett. 96, 180402 (2006). 
2^ Ospelkaus, S., Ospelkaus, C, Wille, O., Succo, M., Ernst, 

P., Sengstock, K., and Bongs, K. Phys. Rev. Lett. 96, 

180403 (2006). 
2^ Catani, J., Barontini, G., Lamporesi, G., Rabatti, F., 

Thalhammer, G., Minardi, F., Stringari, S., and Inguscio, 

M. Phys. Rev. Lett. 103(14), 140401 (2009). 
2^ Bermudez, A., Goldman, N., Kubasiak, A., Lewenstein, 

M., and Martin-Delgado, M. New J. Phys. 12, 033041 

(2010). 
2^ Becker, C, Soltan-Panahi, P., Kronjager, J., Dorscher, S., 

Bongs, K., and Sengstock, K. | arXiv: 091 2. 364 ^1 ? to appear 

in New J. Phys. (2010). 



